Intrinsic intensity fluctuations in random lasers 
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We present a quantitative experimental and theoretical study of intensity fluctuations in the 
emitted light of a random laser that has different realizations of disorder for every pump pulse. A 
model that clarifles these intrinsic fluctuations is developed. We describe the output versus input 
power graphs of the random laser with an effective spontaneous emission factor (/3 factor). 

PACS numbers: 42.25.Dd, 42.55.Zz 



I. INTRODUCTION 

In 1968 Letokhov wrote his pioneering paper jlj] in 
which he predicted that light amplification through stim- 
ulated emission is possible in a random medium with 
gain. A preeminent experimental demonstration of such 
a random laser was published by Lawandy and coau- 
thors in 1994. Typical phenomena of the random laser 
are: a (low) threshold in the power conversion, spectral 
narrowing, and sharp features ("spikes" j_in the emitted 
spectrum for both picosecond 0, M 0. and nanosec- 
ond pump pulses. Due to their low threshold and 
their ease of production, random lasers are expected to 
be used in many applications, such as coding of clothing 
and detection of dangerous materials 9]. The the- 
oretical search for the underlying principles of random 
laser intensifies, focussing on their statistical properties 
and fluctuations. [IS IH El [3 

Recently, a new random-laser phenomenon was de- 
scribed by Anglos and coauthors ^3- They observed 
shot-to-shot intensity fluctuations in the emitted light, 
which were not caused by fluctuations in the pump 
source. In their experiments, the intensity fluctuations 
occur for nanosecond pump pulses, but not for picosec- 
ond pump pulses. The physical understanding of the 
underlying principle of these fluctuations is relevant for 
applications and our perception of random lasers. 

In this paper we present a quantitative experimental 
and theoretical study of the statistics of these shot-to- 
shot intensity fluctuations of random laser systems that 
have different realizations of disorder for every pump 
pulse, and develop a model that clarifles the existence 
of these fluctuations. Our aim is to provide physical 
understanding of the fluctuations. With such a model, 
experimental conditions can be tailored to control the 
fluctuations. 

This paper is organized as follows. In section the 
experimental setup is described, followed in section IIIII 
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by the experimental observations. A model based on the 
number of laser modes in the random laser is presented 
in section Hvl We compare our model with experimental 
observations in section Ivl 



II. APPARATUS AND SAMPLES 

We start by describing the experimental details of our 
samples and the setup. The random laser consists of 
a suspension of Ti02 particles (mean diameter of 180 
nm) in a solution of methanol and Sulforhodamine B (1 
mmol/liter). The suspension is contained in a fused silica 
capillary tube, with internal dimensions 100 x 2 x 0.2 
mm'^. To characterize the mean free path of this sample, 
we performed an enhanced-backscatter cone experiment 
[T^j l and an escape function experiment E3|- We found 
a transport mean free path of 0.46 ± 0.1 /Ltm at A = 633 
nm. 

The samples are excited by a pump pulse at 532 
nm, provided by an optical parametric oscillator (OPO) 
pumped by a Q-switched Nd:YAG laser (Coherent Infln- 
ity 40-100/XPO). The pump pulse has a duration of 3 ns 
and a repetition rate of 50 Hz. The pump light is focused 
with a microscope objective (water- immersed, NA = 1.2) 
onto the sample (focus area = 12 ± 6 ^m^), reaching an 
intensity in the order of 1 mJ/mm^. The light emitted 
by the random laser is collected by the same objective. 
A small part of the pump light is split off from the input 
beam before the objective. This light is used to excite 
a dye solution (Cresyl Violet) which we use as a marker 
for the pump fluence. This dye solution works as a wave- 
length converter for the pump light. The light emitted by 
Cresyl Violet and by the random laser are recorded at the 
same time with a spectrometer and an intensified charged 
coupled device (resolution ^ 0.3 nm spectral width): in 
one frame both the emission of the random laser and the 
pump marker is recorded. The pump light is filtered out 
of the detection path by use of a colored glass filter with 
a transmission of less than 1% at the wavelength of the 
pump laser. 

As a picosecond source we use a stretched pulse at 532 
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nm of a femtosecond optical parametric amplifier (OPA), 
pumped by a femtosecond Ti: Sapphire laser. We could 
not measure the pump pulse duration directly, but cal- 
culated a duration of 15 ± 0.5 ps from the measured 
bandwidth of the OPA and the configuration of the pulse 
stretcher. Measurements of the random laser and Cresyl 
Violet simultaneously were not performed with the pi- 
cosecond pump pulse. 

III. MEASURED FLUCTUATIONS 

We are interested in the intrinsic fluctuations, i.e. the 
fluctuations of the random laser itself which are not the 
result of fluctuations of the pump laser. A typical single 
shot emission spectrum shows two peaks, see the inset of 
Fig.Q the peak of the random laser output is around 594 
nm, and the peak of the marker Cresyl Violet is around 
625 nm. The intrinsic fluctuations of the random laser 
are investigated by comparing the peak height of the light 
emitted by the random laser with that of the marker. We 
have taken many single shot measurements (> 400), and 
from each of these spectra we determine the two peak 
heights. The data is shown in Fig.C] The corresponding 
correlation coefficient is 0.4, indicating that a large part 
of the fluctuations is intrinsic. In contrast, we expect 
stronger correlation between the spontaneous emission 
regime of the random laser and the pump fluence. To 
check this expectation we compute the correlation coef- 
ficient for the spectral radiance at 573 nm, relatively far 
away from the random laser emission peak, and the pump 
marker. For this situation we find a stronger correlation 
coefficient of 0.62, indicating that these fluctuations are 
mainly due to the pump fluctuations. Although the spec- 
tral radiance at 573 nm is still partially influenced by 
stimulated emission, the effect of the spontaneous emis- 
sion can already be seen from the increase of the corre- 
lation coefficient. 

The intrinsic fluctuation coefficient / is defined as the 
ratio between the standard deviation of the shot-to-shot 
intensity of the light emitted by the random laser, A/, 
and the mean value of the intensity of the emitted light 

7. 



To determine the intrinsic fluctuations from the data 
points in Fig. ^ we take a small band of pump fluences so 
that the fluctuations of the pump laser do not influence 
the outcome. We make a histogram of the accompanying 
peak heights of the spectral radiance of the random laser 
and flt this histogram with a gaussian. From this fit we 
obtain the standard deviation and the mean value of the 
peak heights, leading us to the experimental fluctuations 

/nsxxp = (18±3)%, (2) 

where the error margin corresponds to one standard de- 
viation. 
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FIG. 1: The peak of the spectral radiance of the random 
laser is plotted versus the pump marker. The two gray lines 
indicate the mean values of the two peak heights and are used 
to visualize the correlation. The inset shows a typical single 
shot spectrum, from which the counts of the peak heights are 
determined. 

Due to large technical fluctuations of the pump laser, 
no reliable estimate of the intrinsic fluctuations follows 
from the picosecond pumped random laser experiments. 

IV. MODEL 

In this section we present a model for the origin of the 
intrinsic fluctuations of random lasers. The approach we 
take here is based on the concept of pseudo modes l^l- 
Pseudo modes are single frequency eigen modes, solutions 
of Maxwell equations. These modes have an eigenfre- 
quency and decay by leaking to the outside world. This 
leakage is characterized by a decay time Tc. Each mode 
can be a laser mode, depending on the decay time (also 
referred to as dwell time) of that mode and the gain time 
of the sample. The decay time is the time that light 
is inside the sample due to diffusion, and the amplitude 
gain time is defined as the time after which the ampli- 
tude is increased with a factor e. If the decay time of a 
certain pseudo mode is longer than the gain time of the 
system, that pseudo mode is a laser mode. We assume 
the mode volume is equal for each (random) mode, and 
that the gain in the sample is homogeneous. The number 
of lasing modes, Ni, is given by 

Ni = nN, (3) 

with N the number of pseudo modes and tt/ is a random 
variable ranging between and 1. We define 

Pi = , (4) 

with pi the probability for lasing in a pseudo mode, and Wi 
is TT; averaged over realizations of the disorder. The emis- 
sion power of the different lasing modes is equal. The 
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detected emission intensity of the random laser is almost 
entirely due to the lasing modes in the system and is as- 
sumed to be proportional to the number of lasing modes. 
The intrinsic fluctuations / can be determined 
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Ni 
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In the last step of Eq. © , we have assumed a binomial 
distribution of Ni, which results for the limit of N to 
infinity to a gaussian or normal distribution. The stan- 
dard deviation, a, is given by \fWi. We combine Eqs. © 
and 10} , and insert the result in Eq. ||SJ) to obtain 
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In section IIV Al we will show how to calculate pi from 
a fit to experimental data. An elaboration on the calcu- 
lation of N is presented in section llVBI 



A. Determination of the probability of lasing 

When we combine Eqs. Q and Q the probability of 
lasing is given by 



N' 
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The probability of lasing can be calculated via the distri- 
bution of the decay times. The integral of the distribu- 
tion of the decay times P from the gain time to infinity 
will give the probability pi for a certain mode to be a 
laser mode: 



Pi 
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The distribution of the decay times for a 3D diffusive 
medium is not known. We therefore use the distribution 
of the phase delay times 0, which is expected to be 
close to the distribution of the decay times. 
The gain time is given by 



Co 



where n' is the real part of the refractive index, 
amplitude gain length ig is given by 
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where Ue is the stimulated emission cross section of a 
molecule, and pcxc is the density of molecules in the ex- 
cited state in the sample. 

If the pump power is large enough, the gain in the sys- 
tem can be saturated. In the case of saturation the gain 
length will not be decreased any more when the pump 



power is increased. From Eq. (|10() we can find a lower 
bound for the gain length ig i,, and thus an indication 
of gain saturation, when one assumes that all the dye 
molecules in the medium are in the excited state 
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where p is the density of dye molecules in the sample. 

We want to determine the probability of lasing directly 
from our experiments. When we examine Eq. Q, we see 
a similarity between the definition of pi and the sponta- 
neous emission factor of a laser, the (3 factor 19]. The 
singlemode (3 factor, defined as the fraction of sponta- 
neous emission that contributes to lasing is given by |23| 
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This Psra factor appears in the four-level rate equations 
for a singlemode laser [20j 



dNi{t) 

dqjt) 
dt 



(13a) 
(13b) 



with TVi the number of excited molecules in the medium, 
q the number of photons in the lasing mode, Pl the pump 
rate (in photons per second), r the spontaneous emission 
life time of the dye and Tc the cavity decay time. 

In a random laser many random modes contribute to 
the laser oscillation. However, for our consideration only 
the average behavior is relevant. In general, to describe 
a multimode laser one has to write an equation for ev- 
ery mode and couple the different mode equations to the 
equation for the population. Only two small changes to 
the singlemode rate equations are necessary, if the behav- 
ior of the multimode laser can be simulated by the be- 
havior of a singlemode laser. We simply replace /3sm and 
Tc in the rate equations H13|l by the effective parameters 



Pram and Tc 



Since we are interested in the average 



behavior, we will use this simplified approach. We will 
prove that this approach is valid in many situations. 

For the mean value of the cavity decay time, given by 
Tc.mm, we usc the mean value of the distribution of the 
decay times. To determine this distribution we calcu- 
late the solution of the diffusion equation for a slab with 
thickness L, with a source positioned in the middle of 
the slab. From this solution of the diffusion equation the 
electric field correlation is derived and the mean value of 
the phase delay times follows from the Taylor expansion 
of this correlation. 21] The mean value of the phase delay 
time, and thus the mean value of the cavity decay time, 
is given by 



1X2 

8 D 
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The diffusion constant D is given by Co£/{3n'), with £ the 
transport mean free path. 22] For the effective parameter 
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/3mm we take 



/3 

Mmm — J 



(15) 
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where we assumed that ah modes contribute equahy. 
Comparing Eqs. l(7| and ifT^ we see that /3mm is equal to 
Pi- 

In the continuous wave hmit for a singlemode laser in 
steady state a formula can be analytically derived from 
the rate equations (|13() that describes the relation be- 
tween the output and the input power of a laser. 
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/3s, 
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2/5sm 
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If one uses a pulsed pump this should be generalized to 
include time dependence, and an analytic solution is no 
longer available. However, as we will show, for a wide 
range of parameters Eq. (|16|l still describes the threshold 
behavior very well, even for a detector that integrates the 
output power. We will use Eq. (|16|l for the integrated 
power of a multimode laser, with the replacements of the 
parameters as described above. For the use of Eq. Hlt)|) 
in pulsed experiments we replace the parameter /3mm by 
an effective parameter Pcff- 

To examine the applicability of Eq. (|16|l to experi- 
ments with a pulsed pump, we calculate with the rate 
equations (|13|l several output versus input power graphs. 
We use input parameters relevant to our experiment and 
vary the pump pulse duration and /3sm- The pump pulse 
is modelled by a gaussian. To the output versus input 
power graphs we fit Eq. H16|l and use as fit parameter 
/3cS- In Fig. 121 we present the calculated /3cS as function 
of /3mm- For a pulse duration tp of 3000 ps, the calcu- 
lated values of Pcff are identical to the input value of 
/3mm- This correspondence implies that for our system 
nanosecond pump pulses can be treated as a continuous 
wave pump. When the pulse duration is shorter, while 
keeping the other parameters constant, we notice a de- 
viation from this straight line to lower values of /3eff for 
the same values of /3mm- This difference between /3cff and 
/3mm increases for shorter pulses. The origin of the dis- 
similarity between /3mm and /3cff is due to the fact that 
the pump is not a continuous wave, but a pulse with a 
finite duration. 

Besides the pump pulse duration and the spontaneous 
life time, there is a third time scale in the rate equations: 
the mean cavity decay time. If the pump pulse dura- 
tion approaches the value of the cavity decay time Tc^mm 
Eq. ifTCjl is no longer a good fit to the output versus in- 
put power graph. This failure of the fit means that f3cs 
is no longer a parameter that can be used to describe the 
experiment and thus that for that experiment there is no 
/3 factor. In short: 

tp >t:$>Tc CW hmit: /3cff = /3mm, 
T > tp > Tc conversion needed: PcS < /3mm, 
T > Tc > tp simplified model fails. 
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FIG. 2: The theoretical values of the effective /3 factor, /3cff, 
numerically calculated from the rate equations, versus the 
input parameter /3mm for different pump pulse durations. We 
choose the input parameters relevant for our experiments: r 
— 3200 ps and Tc — 0.1 ps. In the legend the pump pulse 
durations are listed. 



We conclude that the rate equations and the threshold 
curve (|16|l can be used for a multimode laser, and that 
these formulas can be used under certain conditions for 
a pumped laser. This allows us to extract pi from the 
threshold curve of a random laser. 



B. Calculation of the number of modes 

The second important parameter to calculate the in- 
trinsic fluctuations is N. The total number of modes in 
the system within the relevant frequency bandwidth Auj, 



can be calculated using the formula 
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m 



N = piu;,V) 



Alo 



Snn^V AA 

A3 a7' 



(17) 
(18) 



with p(uj, V) the density of field modes in the cavity 
volume V, AA the full width at half maximum of the 
emission spectrum, and Ac the central wavelength of the 
emission spectrum. All parameters can be deduced from 
experiments, except the volume of the cavity V. In our 
case, the absorption length is much larger than the trans- 
port mean free path, and we can assume a gain volume 
in the form of hemisphere 



with r the radius of the gain volume. 



(19) 
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FIG. 3: Peak of the measured spectral radiance of the random 
laser versus the pump fluence for a random laser pumped with 
nanosecond pulses. A normal threshold behavior is observed. 
The solid line is a fit to our data, with rc,mm = 0.1 ps, and 
the spontaneous emission life time is 3200 ps. From the fit we 
obtain /3cff = 0.07, which in our model equals the probability 
oflasingp;. [Eq. 11161 1 



The intrinsic fluctuations depend on the number of 
modes, and the probability of lasing. By experimentally 
varying the quantities of the parameters, the intrinsic 
fluctuations of the random laser can be controlled. 



V. RESULTS 

We will now calculate with our model the intrinsic in- 
tensity fluctuations of a random laser system pumped 
with nanosecond or picosecond pulses, and compare these 
fluctuations with experimental observations. 



A. Probability of lasing 

We have measured the peak of the spectral radiance 
of the random laser as a function of the pump fluence, 
as shown in Fig. |2| In case of nanosecond pump pulses 
we have seen that pi = (3cS- From the fit of Eq. to 
the experimental data we find directly the probability of 
lasing 

= 0.07 ± 0.03, (20) 

where the error margin corresponds to one standard de- 
viation. 

For the picosecond pumped random laser we have also 
measured the peak of the spectral radiance of the random 
laser as a function of the pump fluence, see Fig. ^ The 
parameter f3cS has to be converted to /9mm (pi)- In our 
picosecond experiment the pulse duration is 15 ± 0.5 ps 




FIG. 4: Peak of the measured spectral radiance versus pump 
fluence for a random laser pumped with picosecond pulses. 
The solid line is a fit to our data, with the parameters r = 
3200 ps and the mean cavity decay time 0.1 ps. From the fit 
we obtain /3off = 0.03. [Eq. lfT?)jl ] 
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FIG. 5: The calculated /3off versus the input parameter of 
the rate equations /3mm is shown. The graph is produced for 
a pump pulse duration of 15 ps, rc,mm = 0.1 ps and r — 3200 
ps. The solid lines connect the data points. We convert our 
f3cB of 0.03 to a pi of 0.09. 



and a relevant conversion graph for /3mm and PcS is pre- 
sented in Fig. |31 With this graph we convert our values 
of (3cs of 0.03 ± 0.006 to 

Pi p, = 0.09 ±0.015, (21) 

where the error margin corresponds to one standard de- 
viation. 
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B. Number of modes 

The calculation of the number of modes for our ran- 
dom laser regime is given by Eq. 118(1 . In our case we 
have n = 1.4837, Ac = 595 nm. The width of the emis- 
sion spectrum above threshold is 4.3 nm. The volume of 
the gain medium is given by the volume of a hemisphere, 
see Eq. ()19|l . The spatial form of the luminescence com- 
ing from the surface of the random laser sample pumped 
with nanosecond pump pulses was recorded with a charge 
coupled device, while wc filter the pump light. We mea- 
sure a circular spot with a mean radius of 5 ± 0.5 /im. 
The number of modes in the nanosecond pumped situa- 
tion is 

Nns = 746 ± 256. (22) 

We did not record the spatial form of the luminescence 
in case of picosecond pump pulses. However we can spec- 
ulate about the gain volume in the picosecond pumped 
system if we have gain saturation. For our random laser, 
the probability of lasing for the picosecond and for the 
nanosecond pumped case are within each other error mar- 
gins, an indication of gain saturation which wc shall prove 
now. 

We calculate the lower bound of the gain length with 
Eq. Hll|) . The density of molecules in the random laser is 
5.4 X 10^^ molecules m"'^, and the stimulated emission 
cross section is 4 x 10~^° m^, both with a 10% error, 
leading to a lower bound limit of 

Ig^b > 75 ^m. (23) 

For the nanosecond pumped random laser we found a pi 
of 0.07 ± 0.03, see Fig. |21 The gain length at which pi 
equals 0.04 is 

IgAc = 60 fim. (24) 

The value of £g^b,c is lower than the bound limit we calcu- 
lated. Although this does not indicate a serious discrep- 
ancy as the error margins on pi are only one standard 
deviation, the assumption that the phase delay time dis- 
tribution equals the cavity decay time distribution should 
be investigated. Both gain lengths above are almost 
equal, proving that we are near the saturation regime. 
The threshold pump fluence is a factor 100 higher in 
the case of the picosecond pump pulse, and the pulse 
duration of the picosecond pump pulse is a factor 100 
shorter. Since both pump lasers pump the sample to 
gain saturation, the gain volume of the picosecond pump 
pulse is larger than the gain volume of the nanosecond 
pump pulse. The advantage of using two different pump 
lasers on one sample, is the change in the number of laser 
modes. 

The gain volume for the picosecond case in saturation 
is at minimum a cylinder with a radius equal to the lu- 
minescence spot of the nanosecond pumped random laser 
and a length equal to Ld- The length Ld is the length 



that light travels from a point source inside a diffusive 
medium 

Ld = (25) 

with t the pulse duration of the point source. In our 
system Ld is 22 ± 3 fim. This leads to a total number of 
modes in the picosecond pumped random laser of 

A^ps > 2407 ± 885. (26) 

C. Intrinsic intensity fluctuations of a random laser 

The intensity fluctuations derived from our model are 
given by Eq. ||SJ) 

/„, = (14±5)%, (27) 

where the error margin corresponds to one standard devi- 
ation. This number for the fluctuations is in good agree- 
ment with our experimental observations of (18 ± 3)%. 

The intrinsic fluctuations calculated for a picosecond 
pump pulse for our own system is 

/p, < (6.8 ±1.2)%, (28) 

where the error margin corresponds to one standard de- 
viation. We could not verify this result experimentally. 
Anglos and coauthors have performed measurements on 
a random laser pumped with picosecond pulses. [T^ From 
their paper we can calculate the fluctuations of their sys- 
tem with our model. From the flt to the published output 
versus input power graph, and with a conversion from P^s 
to pi, we flnd for their system pi — 0.05. For the gain 
volume we assumed a cylindrical form, as their excitation 
spot is much larger than their mean free path. We obtain 
fps — 0.01%. From their published measured spectra we 
can flnd an estimation for the measured fluctuations. The 
measured fluctuations are (3.7 ± 1.7)%, where the error 
margin corresponds to one standard deviation. These 
fluctuations include the pump fluctuations. These pump 
fluctuations can be in the order of 3% for a typical pi- 
cosecond laser source. Since the expected intrinsic fluctu- 
ations are much smaller than the pump fluctuations, the 
intrinsic fluctuations cannot be measured and our model 
is not inconsistent with their data. 



VI. CONCLUSIONS 

We have developed a model based on quasi modes that 
predicts the fluctuations of the output power of a ran- 
dom laser pumped with either nanosecond or picosec- 
ond pulses. For the system pumped with nanosecond 
pulses we computed fluctuations of (14 ± 5)%. This is 
in good agreement with our experimental fluctuations of 
(18 ± 3)%. For a system pumped with picosecond pulses 
we calculated fps = 0.01%, too small to be observed, for 
the system measured by Anglos and coauthors 14). 
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The difference in intrinsic fluctuations between pi- 
cosecond and nanosecond pumped random lasers is well 
described by our model and the predictions are identi- 
cal to our observations for a nanosecond pumped ran- 
dom laser and published observations for a picosecond 
pumped random laser. Our model can be used to tailor 
experimental conditions in such a way as to control the 
intrinsic fluctuations of a random laser system. 
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